Primordial Star Forming Regions in a CDM 

Universe 

Yu Zhang, Michael L. Norman, Peter Anninos, Tom Abel 

Laboratory for Computational Astrophysics 
University of Illinois at Urbana- Champaign 
405 N. Mathews Ave., Urbana, IL 61801 



Abstract. 

We developed a three-dimensional 2-level hierarchical cosmological code with 
a realistic and robust treatment of multi-species non-equilibrium ionization and 
radiative cooling processes, and use it to investigate primordial star forming 
regions that originate from high-cx perturbations in a standard CDM domi- 
nated cosmology. We find it is possible to produce gravitationally bound and 
cooled structures at very high redshift (z ~ 40) with baryonic masses as small 
as ~ 10 3 Mq. The molecular hydrogen formation in these small scale structures 
follows very well the analytical predictions of Abel (1995) and Tegmark et al. 
(1996). We also discuss the minimum mass that cosmological structures must 
have in order to be able to cool and collapse. 

I INTRODUCTION 

Models for structure formation are based on the growth of small primordial 
density fluctuations by gravitational instability on a homogeneously expand- 
ing background universe. Depending on the nature of the dark matter and 
whether the primordial fluctuations were adiabatic or isothermal, the first 
mass scale to collapse could be as small as one solar mass (very heavy cold 
dark matter particles) or as high as ~ 1O 12 M (hot dark matter scenarios). In 
CDM cosmogonies the fluctuation spectrum at small wavelengths has only a 
logarithmic dependence for mass scales smaller than ~ 1O 8 M , which indicates 
that the small scale fluctuations in this model collapse nearly simultaneously 
in time. This leads to very complex dynamics during the formation of these 
structures, that can be studied accurately only by using realistic numerical 
computations to model the fluid motion and micro-physical processes as well 
as the dark matter component. 

We have recently been able to develop methods that allow us to study the 
problem in three dimensions [5]. We describe this code briefly in section III but 
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first review the process of H 2 formation during small scale structure collapse 
in the following section. A more extensive discussion of our results is given in 
a separate publication [3]. 

II MOLECULAR CHEMISTRY AND COOLING 

The cooling in small scale fluctuations is dominated by the rota- 
tional/vibrational modes of hydrogen molecules. In primordial gas at low 
temperatures ( £ 6000K) molecular hydrogen can not be destroyed efficiently 
unless there is a radiation flux higher than ~ 3 x 10~ 26 erg cm' 2 s~ x in the 
Lyman Werner Bands. Once self-shielding is important, even higher fluxes 
would be needed. The dominating H 2 producing gas phase reaction is the 
dissociative attachment reaction: H~ + H — > H 2 + e~. In the absence 
of an external UV background one can integrate the rate equations to find 
the molecular hydrogen fraction formed during the collapse of primordial gas 
clouds with neutral hydrogen number density nn, temperature T, and initial 
free electron fraction x to be [1,10]: 

hM - / Ha (0) = ^ \n(x n H k r t + 1) = 1(T 8 T L53 ki(t/t° r + 1), (1) 

where hp a, and k r denote the rates for photo- attachment to H and recom- 
bination of hydrogen, respectively. The production of H2 only depends log- 
arithmically on time with a typical time scale of one initial recombination 
time. The temperature dependence is due to the ratio of recombination and 
H formation time scales, which is a measure of the number of electrons avail- 
able to produce H _ . A typical H 2 fraction of ~ 10~ 3 is produced during the 
collapse of structures with virial temperatures greater than 10 3 K. For initial 
(virial) temperatures higher than 6000 K the charge exchange with protons 
will efficiently destroy H 2 , and equation (1) will not be applicable. However, 
during the collapse of clouds with such high virial temperatures the final H 2 
fraction will, nevertheless, be f H2 (T ~ 6000 K) ~ few x 10~ 3 [2]. 

III NUMERICAL RESULTS AND DISCUSSION 

We achieve high spatial and mass resolution with the two-level hierarchical 
three-dimensional code (HERCULES) that we have developed for cosmology 
[4,5]. This code is designed to simulate structure formation in an expanding 
dark matter dominated universe with Newtonian gravity, multi-fluid hydrody- 
namics, radiative cooling, non-equilibrium chemistry and external radiation 
fields. Furthermore, the code independently evolves the following nine species: 
neutral hydrogen H, ionized hydrogen H + , negatively-charged hydrogen H~ , 
hydrogen molecules H 2 , ionized hydrogen molecules H 2 + , neutral helium He, 
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singly-ionized helium He + , doubly-ionized helium He ++ and free electrons 
e~. The 28 most important chemical rate equations (including radiation pro- 
cesses) are solved in non-equilibrium for the abundances of each of the nine 
species. The reaction rates and an extensive discussion of the chemistry model 
are provided in [2]. We have also implemented a comprehensive model for the 
radiative cooling of the gas that includes atomic line excitation, recombina- 
tion, collisional ionization, free-free transitions, molecular line excitations, and 
Compton scattering of the cosmic background radiation (CBR) by electrons. 

We apply our code to high redshift pre-galactic structure formation and 
evolution, investigating specifically the collapse of the first high-a bound ob- 
jects with total masses in the range 10 5 — 1O 9 M . Our model background 
spacetime is a flat (Qq = 1) cold dark matter dominated universe with Hubble 
constant Hq = 50 km s _1 Mpc -1 , baryonic fraction Qb = 0.06, and a hy- 
drogen mass fraction of 76% . The baryonic matter is composed of hydrogen 
and helium in cosmic abundance with a hydrogen mass fraction of 76% and 
ratio of specific heats 7 = 5/3. The initial data for the baryonic and dark 
matter perturbations is the Harrison-Zel'dovich power spectrum modulated 
with a CDM transfer function and normalized to the cluster scale cr 8h -i = 0.7. 
The data is initialized at redshift z = 100 using Bertschinger's [6] constrained 
realization procedure to construct 3 and 4<r fluctuations in cubes of comov- 
ing length 1024kpc, 512kpc, and 128kpc, with total masses of 7.5 x 10 10 M , 
9.3 x 10 9 M Q , and 1.5 x 10 s M , respectively. 

We reproduced the work of Tegmark et al. 1996 with the same cooling 
function [9] we have used in our cosmological hydrodynamics code so that we 
can compare their findings directly to our 3D numerical results. (We note that 
Tegmark et al. used a modified form of the Hollenbach and McKee (1979) H 2 
cooling function.) Figure la is analogous to Fig. 6 of Tegmark et al. 1996 but 
includes the mass evolution in our numerical results for the 4<r perturbations. 
The dotted lines are found by adding up the total mass M 20 o found in cells 
with dark matter overdensities exceeding 200. It is obvious that the use of a 
different cooling function has a very strong influence on the predicted mass 
that can collapse. Although, the quantitative results are very different, the 
shape for these two different M c {z) curves is rather similar at redshifts 30 < 
z < 100. The slopes are consistent with M c oc (l + z)~ 3 / 2 indicating a constant 
virial temperature since T vir oc M 2//3 (l + z). A constant virial temperature in 
turn implies a roughly constant final H 2 fraction given by equation (1). For 
the case in Tegmark et al. , the virial temperature in that regime is ~ 1000 K 
which yields fn 2 ~ 4 x 10 -4 which they argued to be roughly a constant 
universal value which, if exceeded, allows the cloud to collapse at its free fall 
rate. Using the Lepp and Shull H 2 cooling function we find that the virial 
temperature needed to fulfill the Tegmark et al. requirement for collapse is 
m 200 K for redshifts > 30. This translates to a molecular fraction of only 
~ 3 x 10~ 5 . Our simulations, however, show that, although we used the Lepp 
and Shull cooling function, the H 2 fraction at the time when the baryons are 
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FIGURE 1. (a) Collapsed mass vs. redshift. The dark shaded region depicts the 
mass scale for which the virial temperature equals the CMB temperature. Only 
above the light shaded area labeled with LS83 are structures able to collapse. The 
dotted lines show M200 from our 4a 128, 512, and 1024kpc (crosses, squares, tri- 
angles) runs. The dashed upward sloping line is the BBKS CDM spectrum scaled 
appropriately for 4<r peaks. The gray solid downward sloping line is the Jeans Mass 
at 187r 2 the background density. The dot-dashed line shows the original delimiting 
line computed by Tegmark et al. (1996) which was based on a modified form of the 
Hollenbach and McKee (1979) H2 cooling function, (b) The molecular hydrogen 
fraction vs. redshift in the densest zone of the 4<rl28 (open circles) and the 4<t512 
simulations (triangles). The thick solid lines is the solution (1) initialized with the 
appropriate free electron densities, initial H2 fraction, and temperature. 

collapsing into the DM potential wells is also about ~ 5 x 10~ 4 . Hence, we 
find roughly the same critical H2 fraction as derived by Tegmark et al. (1996) 
even with a very different cooling function. 

In Figure lb we test equation (1) against results from the 4<rl28 kpc (open 
circles) and the 4cx512 kpc simulations. The fit is astonishingly good, although 
the initial temperature (or redshift) is somewhat difficult to pick out in this 
case since the heating due to adiabatic compression is slow and the initial 
(virial) temperature is not as well denned as in the case of collapse on larger 
mass scales. Analyzing the time derivative of equation (1) one finds for large 
times that fu 2 ^ T L53 /t which, for the spherical collapse model, translates to 
fn 2 HqM{1 + Zmr) 1 ' 53 when we compare the slopes for different mass scales 
at the present time. This explains that the divergence of the two graphs in 
Figure lb for low redshifts is due to differences in collapse mass and redshift. 

It has been stressed by various authors that early small scale structure might 
influence the entire pregalactic medium and subsequently play an important 
role for structure formation on larger mass scales (e.g. [7]). They are in prin- 
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ciple capable of ionizing a large fraction of the pregalactic medium as well as 
to enrich it with metals. We hope to achieve the required dynamical range in 
future work to estimate the IMF using adaptive mesh refinement techniques 
and so be able to quantify the feedback of collapsing small scale structure. 
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3880 at the NCSA. 
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